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ABSTRACT  This  sets  up  relative  dynamic  models  associated 
with  spacecraft  rendezvous,  does  quantitative  analysis  of 
influences  associated  with  various  types  of  perturbations,  and 
derives  truncation  errors  associated  with  the  processes  of 
setting  up  models.  In  conjunction  with  this,  it  give  calculation 
examples .  Results  are  capable  of  being  used  in  the  analysis  of 
spacecraft  rendezvous  processes  as  well  as  error  corrections. 

KEY  WORDS  Spacecraft  rendezvous  Dynamic  model  Error 
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GRAPHICS  DISCLAIMER 


All  figures »  graphics »  tables  9  equations  9  etc.  merged  into  this 
translation  were  extracted  from  the  best  quality  copy  available . 


INTRODUCTION 


Rendezvous  and  docking  is  one  type  of  high  level  process 
associated  with  astronavigational  activity.  It  is  an  aid  in 
completing  a  considerable  number  of  space  missions — for  example, 
assembly  in  orbit,  increasing  or  exchanging  space  station 
compartment  sections,  experimental  instruments  and  equipment,  as 
well  as  other  useful  loads;  recovering  products,  adding  fuel, 
supplementing  provisions,  materiel,  and  raw  materials;  space 
maintenance  and  rescue  missions,  and  so  on. 

Rendezvous  refers  to  two  or  more  spacecraft  meeting  each 
other  in  orbit  in  accordance  with  predetermined  locations  and 
times.  Docking  refers  to  two  spacecraft  structurally  connecting 
into  one  body  after  iiiceting  each  other  in  orbit.  Speaking  f^om 
the  viewpoint  of  flight  dynamics  and  control,  rendezvous  and 
docking  belongs  to  the  category  of  spacecraft  orbital  control  and 
attitude  control.  In  theoretical  terms  as  well  as  in  terms  of 
technological  realization,  it  is  quite  a  complicated  process  in 
both  cases.  As  a  result,  when  setting  up  dynamic  models  for 
spacecraft  rendezvous,  it  is  not  only  required  that  they  be 
capable  of  comparatively  accurately  scribing  out  the  relative 
motion  processes  of  spacecraft  but  also  that  they  be  convenient 
for  analyses  and  solutions.  Traditional  methods  for  setting  up 
models  often  go  toward  two  types  of  extremes.  One  type  is  making 
use  of  classical  flight  dynamics  principles  to  obtain  simple 
formulae,  yet  it  is  still  difficult  to  describe  spacecraft 
motions.  The  other  type  is  to  consider  the  influences  of 
multiple  kinds  of  factors.  It  is  still  difficult  to  have 
analytical  results.  It  is  only  possible  to  rely  on  numerical 
value  calculations .  This  article  opts  for  the  use  of  traditional 
methods  for  setting  up  models,  quantitatively  analyzing  the 
influences  of  atmospheric  resistance,  solar  light  pressure, 
global  deviations  from  a  true  sphere,  and  the  gravitation  of  sun 
and  moon.  After  that,  it  derives  truncation  errors  associated 
with  the  process  of  setting  up  models,  orbital  errors,  and  so  on, 
convenient  to  estimating  model  accuracies,  offering  a  possibility 
for  precision  positioning. 

1  THE  SETTING  UP  OF  RELATIVE  MOTION  EQUATIONS  FOR  RENDEZVOUS 

Assuming  that  the  target  spacecraft  is  moving  in  the  near 
circular  orbit  I,  the  target  orbital  coordinate  system  OXYZ  is 
adopted  as  the  dynamic  coordinate  system.  Its  origin  point  O  and 
the  center  of  mass  of  the  target  spacecraft  are  firmly  connected. 
In  conjunction  with  this,  it  follows  its  motion  along  orbit  I. 

The  X  axis  is  within  the  orbital  plane.  In  conjunction  with 
this,  it  is  opposite  to  the  direction  of  the  velocity  of  motion. 
The  directions  of  the  Y  axis  and  the  vector  '‘2  associated 
with  the  target  spacecraft  center  of  mass  in  inertial  space  are 
the  same,  as  shown  in  Fig.l. 

In  addition,  assume  that  the  pursuing  spacecraft  center  of 
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mass  vector  in  inertial  space  is 


From  Newtonian  laws 


(3p 


Fig.l  Orbital  Coordinate  System 

Key:  (1)  Target  Spacecraft  (2)  Pursuing  Spacecraft 


as  well  as  dynamics  relationships,  one  gets 
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(1) 


In  the  equation,  i  i 

that  the  following  spacecraft  is  moving 
target  spacecraft,  in  comparison  with 
small  quantity.  Moreover, 


.  Due  to  the  fact 
in  the  vicinity  of  the 


IS  a 


(2) 


As  far  as  equation  (2)  is  concerned, 
the  series,  equation  (1)  is  projected 
OXYZ,  and  one  then  obtains 


expanding  in  accordance  with 
toward  coordinate  system 
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x  —  2(oy  m 

y  +  2o}x  —  3a)  y  ^  T ^  /  m 

V  -¥<0  2=7 ^/m 


(3) 


The  equation  set  in  question  will — within  a  certain  range — 
describe  the  relative  motions  of  target  spacecraft  and  pursuing 
spacecraft. 

2  MODEL  ERROR  ANALYSIS 

(1)  The  Influences  of  Various  Types  of  Perturbations 

During  the  process  of  deriving  equation  (3),  spacecraft  are 

handled  as  problems  where  they  only  act  as  two  bodies  and  have 
not  involved  all  such  influences  as  atmospheric  resistance,  solar 
light  pressure,  the  gravitational  forces  of  the  sun  and  moon,  as 
well  as  global  deviations  from  a  true  sphere,  and  so  on. 

Moreover,  as  far  as  spacecraft  carrying  out  rendezvous  are 
concerned,  they  are  generally  low  orbit  spacecraft.  Densities 
and  velocities  are  both  comparatively  large.  Therefore, 
atmospheric  resistance,  among  perturbation  quantities,  occupies  a 
relatively  large  part.  It  will  consume  most  of  the  spacecraft 
energy. 

When  spacecraft  are  located  outside  the  shadow  of  the  earth, 
they  undergo  the  effects  of  solar  radiation  pressure.  Light 
pressure  direction  is  from  the  sun  pointing  toward  the 
spacecraft.  Earth  shadow  effects  in  light  pressure  perturbations 
are  capable  of  giving  rise  to  abrupt  changes  in  orbital  altitude. 
The  order  of  magnitude  can  reach  several  tens  of  kilometers . 
Besides  this,  the  gravitational  forces  of  the  sun  and  moon  will 
also  produce  a  definite  influence  on  relative  motions  of 
spacecraft.  The  order  of  magnitude  is  103  smaller  than  the 
influence  of  the  shape  of  the  earth. 

(2)  Circular  Orbit  Approximations 

During  the  process  of  derivation,  elliptical  orbits  are 
taken  as  approximations  of  circular  orbits,  ignoring  the 
influence  of  w  and  the  influence  of  rates  of  flattening  in 
series  expansions.  On  the  basis  of  the  maximum  values  of 

and  co’--y  — through  numerical  value  methods — it  is 
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to  calculate  the  influence  of  degrees  of  orbital  non-roundness . 


(3)  Model  Approximation  Treatment  Errors 

Perturbation  errors  achieve  correction  by  the  introduction 
of  perturbation  forces  into  control  quantities.  In  this  section, 
we  discuss  the  errors  which  approximation  treatments  bring  with 
them  in  the  process  of  deriving  formulae. 

((1))  In  the  derivation  process,  take 


and  approximate  it  as 


(r,/r, )  =  1  -~(y/rj 


(4) 


The  truncation  error  Rjp)  can  be  calculated  as  being 

|iJ^(p)|<(3f/2)|p|' 


In  the  equations. 

In  equation  ( 6 ) , 


M  $=  maxi /^(p)  I 


(5) 

(6) 
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Therefore, 


/®(P)  = 


Due  to  y  «  r2 ,  p  «  r2  ,  and  after  omitting  high  order 
terms,  one  has 


« 

r  ^  ^2 
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Finally,  one  obtains 


or 


3  P 
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( ( 2 ) )  Orbital  Errors  Produced  by  Rx  (p) 

Assuming  T=0,  from  equation  (1),  one  has 

^  +2(w  x^)+iox(ax^+a}x  p  =  p - T - 

St  \  St  J  r. 


With  regard  to 


it  is  possible — after  considering  error  Rx{p)  — to  expand  into 


(^) 


1  =1-3—  +rap) 


Equation  (8)  can  be  written  to  become 


+  2  f  0)  X  +  C0  X  (aJ  X  ^  +  0)  X  p 

SA  ^  ^ 


(9) 


Comparing  the  equation  above  and  equation  ( 8 ) ,  it  is 
possible  to  see  that  the  error  below  exists 


As  far  as  conclusions  when  direct  application  is  made  of  existing 
control  forces  are  concerned,  the  orbital  coordinates  are  then 
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In  this,  X,  y,  and  z  are  the  results  when  there  are  no  control 

forces  in  the  first  section. 

Thus,  one  obtains  the  error  formulae  as 


2a  3a. 


Ax  -  -f  -  2-i  sin(a)0 - r « - ^  ^ 

CO  <u  <0 

Ay ^  -  2^sin(£fl<)  -  -|-cos(fflO  -  2  « 


(10) 
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3  CALCULATION  EXAMPLES  AND  CONCLUSIONS 


Use  is  made  of  the  results  derived  in  this  article  to 
analyze  orbital  errors  which  a  certain  spacecraft  is  capable  of 
producing.  Assume  that  the  orbital  parameters  of  the  pursuing 
spacecraft  are:  major  hemiaxis  a  =  6  706.221  km,  eccentricity  e 
=  0.024  578,  orbital  angle  of  inclination  i  =  65",  ascending  node 
celestial  diameter  fi  =  30". 

When  the  relative  distances  of  the  two  spacecraft  are  500km 
and  100km,  truncation  errors  Rx(p)  are,  respectively,  8.34m  and 
0.31m.  In  the  first  type  of  situation,  the  orbital  error  given 
rise  to  is 

Aj:-80  883.84[aj<-sin(coOl;  Ay  =  40  441.92[1 -cos(coO];  Az  =  0 
Due  to  there  existing  in  A;c  a  term  associated  with  a  linear 
increase  as  a  function  of  time,  as  a  result,  divergence  is 
relatively  fast.  Moreover,  Ay  is  then  an  oscillating  term.  Ar 
is  zero  from  beginning  to  ena.  Fig. 2  ~  Fig. 4  are  computer 
simulation  results. 
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Fig. 2  X  Direction  Error  Curve 


Fig. 3  y  Direction  Error  Curve 
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Fig«4  2  Direction  Error  Curve 


From  derived  results  and  simulation  calculations,  it  is 
possible  to  know  that,  following  along  with  the  nearing  of 
relative  distances,  truncation  errors  will  tend  toward  zero. 
However,  due  to  the  oscillating  nature  of  orbital  errors™ 
following  along  with  enlargements  of  time  periods — overall  errors 
will  abruptly  increase,  ,As  a  result,  based  on  dynamic  models 
associated  with  analytical  methods,  it  is  only  possible  to 
describe  spacecraft  movements  within  a  certain  range  and  time 
period.  Carrying  track  determination  precision  a  step  further 
is  normally  completed  through  opting  for  the  use  of  ground 
equipment  remote  control  telemetry  methods.  The  authors  explored 
the  possibility  of  making  use  of  GPS  (global  positioning  system) 
in  combination  with  this  article's  models  to  carry  out  autonomous 
track  determination  [ 3 ] . 
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